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ABSTRACT 

This  paper  presents  an  analysis  of  the  Jindalee  Stage  B  Peak  Detector  and  highlights  a 
number  of  deficiencies.  Solutions  are  proposed  for  these  deficiencies  and  the 
improvements  gained  from  their  implementation  are  measured.  These  solutions 
overcome  clustering  about  the  data  sample  points.  In  addition,  enhancements  are 
made  in  the  areas  of:  adaptive  thresholding  to  maintain  a  constant  false  alarm  rate, 
improved  data  storage  to  accommodate  non-vmiform  peak  densities,  and  integration 
with  the  adaptive  clutter  model  to  associate  peaks  with  the  clutter  regions. 
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Advanced  Jindalee  Tracker: 
Enhanced  Peak  Detector 


Executive  Summary 

This  paper  presents  an  analysis  of  the  Jindalee  Stage  B  Peak  Detector  and  highlights  a 
niunber  of  deficiencies.  Solutions  are  proposed  for  these  deficiencies  and  the 
improvements  gained  from  their  implementation  are  measured.  The  result  of  these 
improvements  is  the  enhanced  peak  detector  that  plays  an  important  role  in  the 
advanced  tracking  system  recently  installed  in  the  Jindalee  Over-the-Horizon  Radar. 

The  primary  deficiency  of  the  Stage  B  Peak  Detector  is  the  way  that  it  clusters  the 
detections  about  the  data  sampling  points.  This  phenomenon  may  cause  moving 
targets  to  exhibit  a  staircase  like  progression  through  the  surveillance  region.  In  the 
radar  environment,  the  automatic  tracking  system  must  cope  with  stationary  targets 
such  as  transponders.  When  the  detector  causes  bias  it  is  difficult  for  the  tracking 
system  to  discriminate  moving  and  stationary  targets.  This  is  particularly  the  case  for 
slow  moving  targets.  The  enhanced  peak  detector  presented  in  this  report  is 
demonstrated  to  overcome  clustering  about  the  data  sample  points. 

The  Stage  B  Peak  Detector  is  also  shown  to  exhibit  a  high  level  of  error  in  estimating 
the  location  of  target  detections.  This  error  is  apparent  from  the  large  scatter  of  the 
location  of  calibration  detections.  The  calibration  signal  is  not  propagated  by  the 
ionosphere  and  should  display  very  little  positional  scatter,  liie  enhanced  peak 
detector  presented  in  this  report  is  shown  to  dramatically  reduce  the  amoimt  of  scatter 
about  the  calibration  location  indicating  a  more  accurate  detection  process. 

The  two  main  developments  in  the  enhanced  peak  detector  to  overcome  these 
deficiencies  are  the  use  of  pre-whitened  data  and  logarithmic  compression.  The 
enhanced  peak  detector  detects  and  locates  peaks  using  the  pre-whitened  data  and 
only  uses  the  whitened  data  for  signal  strength  estimation.  This  reduces  the  location 
error  and  also  reduces  the  incidence  of  peak  splitting  that  occurred  with  the  Stage  B 
peak  detector.  The  result  is  a  more  uniform  distribution  of  noise  peaks  and  reduced 
bias  for  target  peaks. 

The  enhanced  peak  detector  also  has  an  adaptive  threshold  to  maintain  a  constant  false 
alarm  rate.  This  allows  the  peak  detector  to  use  a  lower  threshold  to  obtain  high 
sensitivity  when  clutter  conditions  are  good  without  swamping  the  tracker  with  false 
detections  when  clutter  conditions  are  bad.  To  allow  further  adaptation  to  changing 
clutter  conditions,  the  enhanced  peak  detector  is  integrated  with  an  adaptive  clutter 
model.  This  allows  the  association  of  each  peak  with  a  clutter  zone  and  associated 
statistics. 

To  facilitate  these  new  features  and  allow  for  future  development,  a  new  data  storage 
approach  has  been  adopted.  This  gives  the  enhanced  peak  detector  greater  flexibility 
and  creates  a  simpler  interface  for  the  futiue  development  of  Raid  Size  and  Height 
estimation  algorithms. 


The  impact  of  the  enhanced  peak  detector  is  a  higher  quality  detection  process  to  feed 
the  automatic  tracking  system.  The  enhanced  peak  detector  improves  the  tracker's 
ability  to  resolve  the  ambiguous  velocity  of  radar  targets  as  well  as  providing  more 
accurate  location  estimates  and  smoother  region  hand-over.  By  increasing  the 
uniformity  of  noise  detections  and  associating  clutter  zone  information,  the  enhanced 
peak  detector  helps  the  automatic  tracking  system  significantly  reduce  the  number  of 
false  tracks. 
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1.  Introduction 


Peak  detection  is  used  by  the  Jindalee  OTHR  signal  processing  scheme  to  locate 
target  returns  in  the  azimuth,  range  and  Doppler  (ARD)  data  produced  from 
processing  the  radar's  receiver  output.  The  spatial  resolution  of  an  OTHR  is  quite 
low  (of  the  order  of  tens  of  kilometers)  so  targets  may  be  treated  as  point  reflectors. 
It  is  desirable  to  locate  these  point  reflectors  with  sub-cell  resolution  for  the  purpose 
of  accurate  tracking.  Peak  detection  incorporates  interpolation  to  achieve  this 
improved  resolution  [1]. 

The  radar  data  processing  is  divided  into  two  stages.  The  first  stage  takes  the  time 
series  output  from  the  radar  receiver  array  and  forms  the  time  series  data  from  each 
range  and  azimuth  cell.  The  second  stage  forms  the  frequency  spectrum  for  each 
range-azimuth  cell.  The  second  stage  also  includes  data  conditioning  algorithms  to 
remove  clutter  and  interference.  The  final  step  in  stage  two  involves  peak  detection 
where  local  peaks  above  a  threshold  are  passed  to  an  automatic  tracking  system  that 
is  part  of  third  stage  of  processing. 

The  peak  detector  used  imtil  1997  was  installed  in  September  1983  with  a 
Probabilistic  Data  Association  Filter  (PDAF)  for  tracking  air  targets.  This  system 
was  limited  in  scope  because  the  Jindalee  Stage  B  OTHR  used  array  processors  for  all 
stages  of  signal  processing.  Upgrades  to  the  PDAF  filter  were  also  limited  by  the 
capacity  of  the  array  processors. 

In  1989,  an  experimental  Joint  PDAF  tracker  was  developed  in  readiness  for 
expected  increases  in  processing  capacity.  This  tracker  never  performed  as  well  as 
the  Jindalee  Stage  B  PDAF  tracker.  Some  of  the  reasons  were  found  to  be  related  to 
the  peak  detector.  Preliminary  tests  showed  that  it  sometimes  formed  multiple 
peaks  from  a  single  target  retiun  (called  peak  splitting)  and  peak  locations  were 
biased  towards  the  centre  of  the  radar  processing  cells. 

In  1992  D  Mudge  conducted  research  into  various  methods  for  removing  the  bias  in 
the  peak  detector.  His  report  Detecting  Targets  with  Sampled  Data  was  not  completed 
before  he  left  DSTO.  S  Davey  resumed  this  work  in  1996  and  this  report  contains  the 
results  on  bias  reduction  by  D  Mudge  with  new  developments  by  S  Davey. 

The  report  starts  with  a  description  of  the  original  peak  detector  and  its  defects.  The 
sections  that  follow  give  the  solution  to  each  defect  with  an  analysis  where 
appropriate.  The  final  sections  give  additional  improvements  to  the  peak  detector 
with  Ihe  results  from  tests.  In  this  paper  the  original  peak  detector  is  referred  to  as 
either  the  Stage  B  or  the  Jindalee  Stage  B  peak  detector.  The  term  IRSU  (Number  1 
Radar  Surveillance  Unit)  is  used  to  refer  to  the  current  radar  system. 
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2.  The  Jindalee  Stage  B  Peak  Detector. 

The  PDAF  tracking  system  uses  peak  detections  to  initiate  and  update  tracks.  The 
peak  detections  are  found  by  firstly  finding  local  maxima,  referred  to  as  peak 
locations,  in  the  ARD  data  after  data  whitening.  The  ARD  data  are  a  three 
dimensional  array  of  Signal-to-Noise  Ratio  (SNR)  values  in  units  of  power  for  each 
range,  azimuth  and  Doppler  cell.  Peak  locations  in  this  array  are  found  by 
comparing  each  sample  point  to  its  immediate  neighbours  in  each  of  the  three 
dimensions.  The  sample  is  declared  to  be  a  maximrun  if  the  SNR  received  at  that 
position  is  greater  than  each  of  its  neighbours.  Referring  to  Figure  1,  the  sample 
labeled  A  is  a  local  maximum  if  it  is  greater  than  all  of  the  points  marked  B. 


Figure  1  ARD  Values  Used  for  Peak  Location  by  the  Stage  B  Peak  Detector 


This  approach  to  peak  location  is  valid  at  all  data  points  except  for  those  along  the 
boundary  of  the  ARD  data.  Here  there  is  no  outer  neighbour  in  azimuth  and,  or 
range,  i.e.  no  shadow  beams  and  ranges.  Therefore  the  approach  illustrated  in 
Figure  1  is  modified  at  the  ARD  data  boimdaries  to  exclude  the  outer  non-existent 
cell  for  peak  location.  In  Doppler  this  problem  does  not  exist  as  the  land  clutter  is  at 
0  Hz.  Here  the  peak  detector  starts  and  ends  4  cells  inside  the  Doppler  cell  range 
that  corresponds  to  doppler  frequencies  from  OHz  to  Waveform  Repetition 
Frequency  (WRF). 

In  the  peak  selection  algorithm,  all  of  the  local  maxima  for  a  particular  range- 
azimuth  cell  are  found  and  sorted  by  sample  SNR  value.  From  these  maxima, 
potential  peaks  are  chosen  by  selecting  the  N  largest  maxima  in  the  range-azimuth 
cell.  Storage  for  a  fixed  number  of  peaks  gives  a  fixed  sized  array  that  allows  the 
peaks  to  be  stored  according  to  their  location  in  azimuth-range  space.  This  allows 
detection  and  tracking  algorithms  to  efficiently  utilise  the  Jindalee  Stage  B  array 
processor.  However  this  limit  degrades  detection  performance  when  there  are  either 
multiple  targets  or  a  target  with  spread  clutter  in  a  single  range-azimuth  cell.  In  a 
test  using  a  nighttime  data  set  that  contained  targets,  this  limit  was  exceeded  in  57 
percent  of  dwells.  If  we  increase  the  munber  of  peaks  in  each  cell  to  overcome  this 
limitation,  a  large  percentage  of  the  peak  array  will  be  empty.  Also  there  is  an 
increase  in  the  processing  time  being  used  to  process  non-data  signals. 
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Once  N  maxima  have  been  chosen,  the  position  of  each  peak  is  determined  by 
interpolating  about  the  position  of  the  local  maximum  using  a  parabolic  curve  fit,  as 
shown  in  Figme  2.  This  reduces  the  peak  detector's  measurement  error  for  target 
tracking.  The  SNR  of  each  peak  is  taken  from  the  sample  value  at  the  local 
maximum.  After  peak  interpolation,  each  peak  is  compared  to  a  fixed  threshold  and 
canceled  if  it  does  not  exceed  it.  The  thresholded  peaks  are  loaded  into  the  peak 
array  and  any  imused  space  is  filled  with  dummy  peaks. 


interpolated 
peak  position 


i- 1  i  i  + 1 


Figure  2  Peak  Interpolation  Using  a  Parabolic  Fit 


> 


As  the  interpolation  process  uses  neighbours  in  all  directions  about  a  data  peak, 
artificial  neighbors  are  generated  for  data  peaks  along  the  edges  of  a  range-azimuth 
boimdary  of  the  data.  The  artificial  neighbours  are  formed  by  using  a  scaled  down 
value  of  the  edge  cell  peak  value. 

The  process  is  summarised  in  the  flow  diagram  shown  in  Figure  3. 
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Figure  3  Jindalee  Stage  B  peak  detection  flow  diagram 


2.1  Stage  B  Algorithm  Defects 

The  defects  of  the  Stage  B  peak  detector  are  best  illustrated  by  the  scatter  plot  in 
Figure  4.  The  scatter  plot  shows  the  distribution  of  interpolated  peaks  in  range  and 
azimuth  from  a  single  radar  region  accumulated  over  188  radar  dwells.  The  data 
contains  a  single  target  and  a  stationary  calibrate  signal.  The  interference  is 
primarily  external.  Hence  a  uniform  distribution  is  expected.  Note  the  tendency  for 
the  peaks  to  be  clustered  around  the  range-azimuth  cell  positions.  A  similar 
clustering  is  found  aroxmd  the  Doppler  cell  positions. 
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The  scatter  plot  highlights  the  stationary  calibrate  signal  midway  between  the 
second  and  third  last  cells  on  the  horizontal  axis.  Instead  of  there  being  peaks 
located  at  a  single  point,  the  peaks  are  scattered  about  this  position.  Analysis  of  the 
number  of  peaks  in  the  region  of  the  calibrate  signal  did  not  show  any  evidence  of 
the  peak  splitting  mentioned  in  the  introduction.  Section  4  contains  an  analysis  of 
peak  splitting  for  a  data  set  containing  multi-mode  targets. 


From  the  above  results  and  as  observed  in  the  introduction,  the  Stage  B  peak 
detector  defects  are  in  the  areas  of : 


•  Peak  interpolation  bias  errors  -  including  edge  effects, 

•  Peak  splitting, 

•  Inflexible  data  structures. 

The  following  sections  deal  with  the  solution  to  these  problems. 
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3.  Peak  Interpolation  Bias 

The  scatter  plot  in  Figure  4  demonstrates  the  biasing  present  in  the  Stage  B  peak 
detector.  The  bias  clusters  peaks  about  processing  cell  positions  and  makes  it  more 
difficult  for  the  tracker  to  correctly  resolve  stationary  and  slow  moving  targets. 

The  main  causes  for  this  problem  are: 

Interpolation  Function  Errors  will  occur  if  the  shape  of  target  returns  does  not 

closely  resemble  the  parabola  used  for  interpolation. 
When  the  target  echo  is  a  point  source,  it  is  as  an 
impulse  in  the  frequency  spectrum.  This  impulse  is 
then  spread  by  the  processing  windows  in  each 
dimension.  Other  factors  that  spread  the  target  rehuns 
are  non-linear  processing  steps,  additive  noise  and 
distortions  in  the  propagation  path. 

Distortions  from  Whitening  Peak  detection  and  interpolation  is  applied  to  data 

after  the  data  whitening.  Data  whitening  scales  the 
ARD  data  to  the  local  noise  floor  so  that  whitening 
data  is  in  SNR.  Because  whitening  is  a  non-uniform 
process  and  because  the  background  noise  is  non- 
uniform,  it  will  distort  the  target  return  and  hence 
increase  the  interpolation  process  error. 

A  number  of  other  possible  causes  for  the  phenomenon  were  investigated  and  these 
are  presented  in  the  following  subsections. 

3.1  Method  of  Analysing  Interpolation  Errors  and  Bias 

The  original  investigation  presented  the  peak  detector  bias  by  forming  a  histogram 
of  interpolated  peak  locations  relative  to  the  processing  cells.  The  histogram  uses 
100  bins  over  the  range  +  0.5  cell.  For  interpolation  with  no  bias,  the  histogram 
should  approximate  a  uniform  distribution.  The  variance  of  the  observed 
distribution  from  a  uniform  distribution  is  used  as  a  figure  of  merit.  The  variance  is 
normalised  for  a  histogram  of  1000  data  points.  The  steps  to  evaluate  the  figxire  of 
merit  are  shown  in  Figure  5.  Here  M  represents  the  actual  number  of  data  points 
used  for  the  histogram  and  N  is  the  niunber  of  bins. 
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Figure  5  Normalised  variance  calculation 

3.2  Bias  Correction  by  Data  Compression 

Because  a  parabola  gives  a  simple  interpolation  algorithm,  the  approach  taken  to 
investigate  the  centering  was  to  alter  the  data  values  so  that  the  shape  of  returns 
gives  a  better  match  to  a  parabola.  The  best  results  came  from  compressing  the  data. 
This  action  maintains  the  peak  location.  The  compression  techniques  investigated 
were: 

a)  Square  root  the  data  array  values  before  interpolation.  Because  this  reduced 
the  bias,  tests  were  also  made  with  4*  and  8*  roots. 

b)  Conversion  of  the  data  values  to  decibels  (i.e.  SNR  values)  before 
interpolation. 

Data  compression  was  tested  on  both  simulated  data  and  actual  OTHR  data.  The 
results  set  out  in  Table  1  and  Table  2  show  the  histogram  variances  for  each  case 
respectively.  It  can  be  seen  that  both  the  8*  root  and  the  decibel  techniques  provide 
significant  improvement  in  performance.  The  best  performance  for  OTHR  data  is 
achieved  with  the  decibel  technique. 


8*^  root 

4*^  root 

Square  root 

dB 

Baseline 

Azimuth 

1091 

1170 

1298 

916 

2104 

Doppler 

999 

989 

998 

965 

1269 

Range 

917 

967 

1228 

883 

1562 

Table  1  Histogram  error  variances  for  data  compression  of  simulated  data 
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8‘*'  root 

root 

Square  root 

dB 

Baseline 

Azimuth 

2270 

2668 

3467 

1987 

5256 

Doppler 

229 

279 

401 

189 

719 

Range 

956 

1153 

1625 

770 

2744 

Table!  Histogram  error  variances  for  data  compression  of  OTHR  data 

While  the  8*  root  and  decibel  techniques  provide  similar  results,  the  8*  root 
technique  is  much  slower.  This  is  due  to  the  implementation  of  exponents  by 
compilers.  The  algorithm  for  taking  the  n*  root  of  a  number  is  shown  below. 


y  =  exp(-n  x  log(x))  =  y  =  Vx 


The  algorithm  contains  an  exponential  operation  and  a  logarithm  both  of  which  are 
quite  computationally  expensive.  The  exponential  and  logarithm  are  approximately 
the  same  speed  for  operation  and  so  the  decibel  technique  takes  approximately  half 
the  time  of  the  8*  root.  This  makes  the  decibel  technique  much  more  desirable. 

3.3  Bias  Correction  by  Look-up  Table 

Another  method  for  correcting  the  bias  is  to  use  a  look-up  table  to  adjust  for  systematic 
errors  in  the  interpolation  process.  The  interpolation  error  is  a  function  of  SNR  with 
the  largest  errors  at  low  SNR.  To  test  the  viability  of  the  look-up  table  method,  the 
table  values  were  chosen  to  correct  the  interpolated  position  of  a  lOOdB  target. 


Look-Up 

Look-Up  on  8**^  Root 
Data 

Look-Up  on  dB 
Data 

Baseline 

Azimuth 

1979 

2494 

2599 

5256 

Doppler 

176 

207 

223 

719 

Range 

925 

1114 

1154 

2744 

Tables  Histogram  variances  for  look-up  correction  of  OTHR  data 


Table  3  gives  the  variances  using  corrections  from  a  look-up  table.  The  performance  is 
similar  to  data  compression  using  the  log  function.  Also  listed  in  Table  3  are  the 
variances  when  data  compression  with  the  8*  root  and  log  functions  is  used  in 
conjimction  with  a  table  look-up.  For  this  case  there  is  a  sUght  degradation  in 
performance. 

The  look-up  table  method  requires  fewer  operations  than  data  compression  because  it 
only  involves  a  table  interpolation  and  subtraction  per  interpolation  dimension. 
Improvement  in  performance  is  expected  by  adding  SNR  to  the  table  but  this  also  adds 
more  operations.  One  feature  of  the  look-up  table  approach  is  the  need  to  have  tables 
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for  each  significant  change  in  signal  processing  weights.  A  related  issue  is  the 
unknown  sensitivity  of  the  look-up  table  to  signal  distortions  introduced  by  the 
ionospheric  medium.  Because  of  these  issues  and  because  the  simplified 
implementation  gave  similar  performance  to  data  compression  with  a  log  function,  the 
look-up  table  is  not  developed  any  further. 

3.4  Bias  Introduced  by  Peak  Thresholding 

The  Stage  B  peak  detector  applies  thresholding  to  the  data  value  at  the  sample  point, 
i.e.  the  radar  processing  cells.  Consider  a  target  peak  with  an  SNR  value  just  above 
the  detection  threshold.  If  the  position  of  the  signal  maximum  is  midway  between 
two  sample  points,  the  sample  values  will  be  lower  than  the  threshold. 
Alternatively,  if  the  maximum  is  close  to  a  sample  point,  the  sample  will  exceed  the 
threshold.  This  effect  is  sometimes  called  scalloping  loss.  Hence  signals  just  above 
the  threshold  and  occurring  near  the  sample  points  are  more  likely  to  be  passed  by 
the  peak  detection  routine  than  similar  signals  near  the  midpoint  between  samples. 
This  is  illustrated  in  Figtue  6.  The  left  return  will  give  a  peak,  but  the  right  one  will 
not. 


The  way  to  remove  this  possibility  is  to  interpolate  in  SNR  before  applying  the 
threshold  for  accepting  peaks.  Interpolation  has  the  added  advantage  of  removing 
SNR  fluctuations  introduced  as  the  target  moves  between  the  radar  processing  cells. 

SNR  interpolation  is  added  by  calculating  the  maximum  of  the  parabola  fitted  to  the 
sample  point  and  its  neighbours.  The  threshold  can  then  be  applied  to  the 
interpolated  peak  value.  In  this  way,  the  same  thresholding  criteria  is  applied  to  all 
peaks  irrespective  of  their  position  relative  to  the  sampling  points,  and  a  more  even 
distribution  should  result  for  signals  near  the  detection  threshold. 

3.4.1  Multidimensional  Interpolation 

We  need  to  extend  the  above  process  of  interpolation  for  quadratics  of  a  single 
variable.  We  want  to  interpolate  the  amplitude  as  a  function  of  three  variables: 
range,  azimuth  and  doppler  shift.  It  is  fairly  simple  to  do  this  if  we  model  the  target 
return  as  shown  in  equation  (1). 
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A(r,a,d)  -  F(r)  xG(a)xH(d) 


(1) 


That  is,  we  assume  that  the  amplitude  function  of  the  target  return  is  separable. 
With  the  model  in  this  form,  we  can  find  the  peak  in  each  dimension  independently 
(this  is  the  process  used  for  position  interpolation).  We  may  calculate  the  value  of 
the  maximum  (the  interpolated  peak)  using  equation  (2). 


A(ro,ao,do)  x 


^  ^  HidJ 

Fir^)  G(flo)  H(do). 


(2) 


where  (r^,  a^,  dj  are  the  coordinates  of  the  maximum  and  a^,  d^  are  the 
coordinates  of  the  ARD  data  peak.  The  expression  in  brackets  is  a  correction  factor 
applied  to  the  sample  value  to  obtain  the  local  maximum.  We  assume  that  it  is  valid 
to  approximate  the  unknown  functions  F,  G,  H  near  the  sample  maximum  using 
parabolic  interpolation. 

Each  of  the  ratios  to  the  right  of  equation  (2)  are  calculated  separately.  The  first 
ratio,  F  (r„  )/F  (tq  )  is  written  as 


F(r„)  ^  F{rJxG(.ao)x  H(d,)  ^  Air„,a„d,)  (3) 

F(rf^)  f  (ro)  X  G(flo)  X  H(do)  A(ro,ao,do) 

where  the  term  A(r^,  Ug,  dg)  is  the  peak  value  obtained  by  applying  one  dimensional 
parabolic  interpolation  in  range.  Similar  expressions  are  used  for  the  azimuth  and 
doppler  ratios  in  equation  (2).  The  ratios  are  then  substituted  into  equation  (2)  to 
obtain  the  interpolated  peak  value. 


n  ^  ^  ^  A(rQ,flQ,d^) 

m  r'*m  m  '  ~  - - Z — Tj - 

Mro,ao,do) 


(4) 


Note  that  the  form  of  the  interpolation  equation  is  the  product  of  one  dimensional 
interpolations.  This  is  due  to  our  initial  assumption  that  the  target  return  provides  a 
separable  function. 

3.4.2  Effect  of  Amplitude  Interpolation 

The  amplitude  interpolation  technique  was  included  with  the  Stage  B  peak  detector 
and  the  effect  that  it  had  on  recorded  OTHR  data  was  tested.  There  was  a  significant 
rise  (17%)  in  the  number  of  peaks  with  a  fixed  threshold.  However,  the  additional 
peaks  were  still  clustered  about  the  sample  points.  Only  a  small  improvement  in  the 
distribution  was  observed.  This  is  illustrated  by  the  variances  calculated  from 
recorded  data  in  Table  6. 
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Interpolated  Baseline 

Azimuth 

Doppler 

Range 

4452  5256 

520  719 

2159  2744 

Table  4  Histogram  variances  for  amplitude  interpolation  on  OTHR  data. 


The  process  of  interpolating  the  peak  SNR  moves  the  peak  probability  density 
function.  It  increases  the  mean  and  modifies  the  variance  of  the  distribution.  The 
distribution  of  interpolation  corrections  is  approximately  a  gamma  distribution  as 
shown  below  in  Figure  7. 


Figure  7  Interpolation  Correction  Distribution 


3.4.3  Implementing  a  Dynamic  Peak  Threshold 

Amplitude  interpolation,  the  number  of  targets  and  OTHR  clutter  variations  change 
the  munber  of  peaks  above  a  fixed  threshold.  To  achieve  a  constant  false  alarm  rate, 
a  dynamic  threshold  is  applied  to  the  peaks.  The  method  chosen  to  implement  the 
dynamic  threshold  involves  forming  a  cumtilative  distribution  fimction  of  the  peak 
SNR  values.  This  distribution  has  a  resolution  of  say  0.5  dB  and  the  threshold  is 
chosen  by  selecting  the  SNR  with  the  required  number  of  peaks. 
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The  cumtilative  distribution  requires  a  separate  array  of  peaks  (that  will  not  be  used 
in  data  output)  with  all  of  the  unthresholded  peaks  in  the  dwell.  As  a  peak  is  found, 
the  cumulative  distribution  array  is  updated.  This  array  is  indexed  by  SNR  and  each 
element  contains  the  number  of  peaks  with  an  SNR  greater  than  the  SNR 
corresponding  to  the  index  value.  For  example,  we  might  form  a  distribution  with 
thirty  entries  in  steps  of  0.5  dB.  A  fabricated  distribution  of  this  form  is  shown  in 
Table  5.  If  we  are  choosing  say  200  peaks  for  each  dwell,  the  threshold  in  this  case 
would  be  chosen  to  be  6.0  dB. 


Table  5  Example  Peak  Distribution 

To  successfully  implement  a  dynamic  threshold,  we  cannot  use  a  peak  array  data 
structure  that  has  a  fixed  number  of  peaks  per  range-azimuth  cell.  A  new  data 
structure  for  use  with  a  dynamic  threshold  is  described  in  Section  5. 

3.5  Bias  from  Peak  Interpolation  on  the  ARD  data  boundary 

The  artificial  neighbours  generated  for  non-existent  data  points  along  the  range- 
azimuth  botmdary  of  the  data  introduce  a  peak  interpolation  bias.  This  increases  the 
measurement  error  of  target  peaks  by  adding  a  bias  that  causes  tracks  to  kink  as  they 
cross-region  boundaries. 

Also,  because  the  peak  detection  process  uses  fewer  neighbouring  points  to  locate  a 
peak,  ie  5  instead  of  6,  the  edge  cells  have  a  higher  number  of  peaks.  This 
contributes  to  an  increase  in  the  number  of  false  tracks  along  the  region  boimdaries. 

The  correction  of  this  problem  involves  the  addition  of  shadow  beams  and  ranges  to 
the  stage  one  and  two  processing.  The  peak  detector  is  then  able  to  use  the  same 
peak  location  and  interpolation  algorithm  over  all  ARD  cells. 

3.6  Bias  from  Sub-Array  Beam  Pattern 

The  receiving  array  beamforming  is  performed  in  two  steps.  The  first  consists  of 
multiple  delay  line  beamformers  that  are  connected  to  form  overlapped  subarrays. 
The  output  from  the  sub  array  beamformers  is  input  to  the  radar  receivers  which 
contain  A/D  converters.  The  binary  data  is  then  used  as  input  for  the  final  step  of 
the  beamforming  process  that  forms  the  simultaneous  fingerbeams.  Figure  8 
illustrates  the  subarray  beam  (dashed  line)  and  the  6  central  fingerbeams  with  2 
shadow  beams.  This  figure  also  shows  the  sidelobes  formed  by  this  two  step 
beamforming  process. 
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Figure  8  Subarray  antenna  pattern 


The  main  feature  of  Figure  8  that  contributes  to  a  bias  in  the  interpolation  is  the 
tapering  of  the  fingerbeam  gain.  This  leads  to  an  azimuth  bias  towards  the  central 
fingerbeam.  The  method  for  correcting  this  bias  is  found  from  the  radar  equation. 
For  the  case  of  an  externally  noise  limited  system  with  a  target  at  a  fixed  location,  the 
radar  equation  gives  the  SI^  proportional  to  the  following  system  parameters. 


SNR  oc 


P,G,gr 

iVo 


Where  P,  is  the  transmitter  power,  G,  is  the  transmitter  gain,  g,  is  the  receiving  array 
gain  and  is  the  noise  power  spectral  density.  For  a  particular  target,  only  the 
receiving  array  gain  varies  because  of  the  fingerbeam  overlap  -  the  transmitting 
power  and  gain  serve  only  as  a  target  illuminator.  So  the  transmitting  array  beam 
shape  does  not  contribute  to  any  bias  in  target  location,  it  only  alters  the  maximum 
target  SNR.  Also  the  noise  power  does  not  depend  on  antenna  gain.  Therefore  the 
only  correction  required  for  the  change  in  gain  is  that  from  the  subarray  beam 
shading. 

This  correction  applies  to  tire  peak  selection  and  interpolation.  Note  that  correcting 
for  the  antenna  gain  alters  the  peak  location  in  azimuth.  The  correction  is  applied  by 
scaling  the  pre-whitening  data  using  a  fixed  array.  The  values  in  this  array  are 
determined  by  the  sub  array  shading  and  are  almost  constant  as  the  steer  direction  is 
changed.  The  whole  of  the  pre-whitening  data  array  is  scaled  before  peak  detection. 

3.7  Recommendations  for  Correction  of  Interpolation  Bias 

From  the  above  analysis,  reduction  of  the  interpolation  bias  errors  in  the  Stage  B 
peak  detector  requires  the  introduction  of  the  following  changes. 
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a)  Data  compression  with  a  log  fimction  to  give  SNR  in  dB. 

b)  Amplitude  interpolation.  (Note  that  amplitude  interpolation  does  not 
significantly  reduce  bias,  but  it  does  improve  the  accuracy  in  the  measured  SNR.) 

c)  Addition  of  shadow  beams  and  ranges  to  ARD  data. 

d)  Correction  for  the  fingerbeam  gain  taper  by  the  subarray.  (Note  that  this  change 
does  not  contribute  to  the  bias  observed  in  noise  peaks  because  the  noise  is 
isotropic.  However,  it  adds  an  error  to  the  interpolated  azimuthal  position  for 
true  target  returns.) 

With  these  changes  implemented,  the  new  variances  for  tests  performed  on  OTHR 
data  are  given  below  in  Table  6.  The  relative  improvements  are  different  for  each 
dimension.  This  is  probably  due  to  the  fact  that  different  windows  are  applied  for 
each  dimension  in  the  signal  processing. 


Proposed 

improvements 

baseline 

Azimuth 

1021 

5256 

Doppler 

54 

719 

Range 

228 

2744 

Table  6  Histogram  variances  with  recommended  improvements  implemented 
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4.  Peak  Splitting 

The  tendency  of  the  Stage  B  peak  detector  to  produce  multiple  peaks  about  a  single 
target  return  is  particularly  evident  for  high  SNR  returns.  The  Jindalee  PDAF 
tracking  filter  assumes  that  a  target  gives  a  single  return.  The  new  multiple  model 
track  initiation  filter  [2]  is  particularly  sensitive  to  multiple  closely  spaced  peaks 
because  they  can  make  velocity  ambiguity  resolution  more  difficult. 

One  cause  for  the  splitting  effect  is  skewness  of  peaks  in  the  ARD  data.  A  mild  skew 
exists  because  of  the  range-Doppler  coupling  ^at  is  introduced  by  processing  the 
waveform  in  two  stages,  i.e.  range  followed  by  Doppler  processing.  Also  skew  can 
be  introduced  by  noise.  Figure  9  demonstrates  how  skewed  data  gives  extra  peak 
detections.  Point  A  is  the  true  peak,  however,  points  B  and  C  will  be  declared  as 
peaks  because  of  their  neighbours.  The  neighbours  used  for  maxima  definition  are 
shown  as  imfilled  circles. 


Another  cause  for  splitting  is  a  result  of  the  peak  selection  being  after  the  data 
whitening  process.  The  Stage  B  whitener  scales  the  ARD  data  by  dividing  it  into 
Doppler  bands  followed  by  range  bands.  In  each  band  a  trimmed  mean  is  calculated 
for  scaling  the  data.  Thus  the  whitener  scales  the  data  in  each  band  by  a  slightly 
different  level,  which  can  lead  to  a  single  broad  peak  being  split  into  a  cluster  of 
peaks  as  shown  in  the  one  dimensional  example  in  Figure  10.  A  similar  splitting  of 
peaks  may  arise  from  multiple-mode  propagation.  In  this  case  the  target  returns 
from  different  propagation  paths  are  displaced  in  Doppler  as  well  as  range.  The 
multi-mode  returns  may  bias  the  whitener  leading  to  distortions  in  the  shape  of  the 
peak  which  will  give  clustering  as  well  as  interpolation  errors. 
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The  solution  to  this  defect  involves  altering  the  method  of  peak  selection  and  its 
relation  to  the  whitening  process.  One  approach  investigated  increased  the  power  of 
each  sample  by  a  weighted  sum  of  its  neighbours.  The  weights  were  based  on  the 
ARD  processing  windows.  After  performing  the  weighted  sum  for  all  ARD  cells,  the 
previous  peak  selection  criterion  was  used  except  that  in  this  case  the  power  array  was 
used  instead  of  the  data  array.  This  approach  would  tend  to  have  a  smoothing  effect 
and  reduce  the  occurance  of  peak  splitting.  The  performance  was  sensitive  to  the 
values  in  the  array  of  weights  and  in  some  cases  the  peaks  were  not  in  the  correct 
location.  So  this  approach  was  not  pursued  any  further. 

The  solutions  given  in  the  following  sections  do  not  alter  the  ARD  data  values  but 
deal  with  alterations  to  the  peak  selection  procedure  and  its  point  of  application  in 
the  processing  sequence. 

4.1  Peak  Location  with  all  Neighbours 

To  stop  skewed  peaks  in  the  ARD  data  forming  multiple  peaks,  the  sample 
neighbourhood  is  extended  to  include  cells  on  the  diagonal  as  shown  in  Figure  11. 
These  extra  cells  increase  the  number  of  cells  tested  for  peak  selection  from  7  to  27. 
This  prevents  the  formation  of  multiple  peaks  as  shown  in  Figure  9. 
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Figure  11  Modfied  neighbourhood  for  local  maxima 


The  peaks  produced  using  the  modified  neighbourhood  have  been  compared  with 
those  using  the  standard  six  neighbour  system.  The  comparison  used  the  same  data 
set  that  produced  the  scatter  plot  in  Figure  4.  This  data  contains  a  single  target 
propagated  via  a  single  mode.  In  this  case  the  number  of  peaks  per  dwell  was 
coimted  for  both  methods  of  peak  location  with  a  fixed  threshold.  Figure  12  shows 
the  percentage  reduction  in  the  total  number  of  peaks  per  dwell  as  a  function  of 
time.  On  average,  the  number  of  peaks  is  reduced  by  at  least  15  percent.  This  is 
clear  evidence  of  a  reduction  in  splitting. 


Figure  12  Reduced  number  of  peaks  from  using  all  neighbours 
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4.2  Peak  Detection  on  Pre-Whitening  Data 

As  noted  previously,  it  is  highly  desirable  to  use  data  before  the  whitening  process 
for  peak  location  and  interpolation.  This  involves  saving  the  pre-whitening  data  for 
later  use  by  the  peak  detector.  The  tracking  system  has  been  foimd  to  work  best 
with  the  SIW  of  peaks  rather  than  the  absolute  power.  For  this  reason  peak  SNR  is 
retained. 

The  data  prior  to  the  whitening  process  is  not  in  a  form  that  is  suitable  for 
determining  the  SNR  of  peaks  because  the  backgroxmd  noise  is  not  white.  It  is 
possible  to  convert  the  data  to  an  absolute  scale  (such  as  dBm)  but  the  thresholding 
for  peak  selection  would  need  to  vary  with  conditions  (particularly  range).  Instead 
the  data  is  normahsed  by  using  the  same  value  to  scale  every  point  to  get  an  average 
noise  floor  of  unity.  Because  the  data  may  not  be  spatially  uiuform,  the  sample 
levels  derived  at  this  point  will  not  give  a  true  measure  of  the  peak  SNR.  The  SNR 
must  be  derived  from  an  approximation  to  the  local  noise  level.  This  is  what  is 
performed  in  data  whitening. 

This  means  that  we  need  the  whitening  data  to  obtain  peak  SNR  values  and  the 
pre-whitening  data  to  locate  and  interpolate  maxima.  This  leads  to  the  flow  diagram 
shown  in  Figure  13. 


PEAKS 


Figure  13  Data  Flow  Diagram  of  peak  detection  on  pre-whitening  data 
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In  Figure  13  the  "locate  and  interpolate  peaks"  uses  the  pre-whitening  data  to  search 
for  peaks  in  the  area  defined  by  the  threshold  crossings.  This  is  more 
computationally  efficient  than  searching  all  of  the  pre-whitening  data  for  peaks  and 
then  removing  those  peaks  that  did  not  pass  the  threshold  test.  Here  the  peak 
location  uses  the  all  neighbour  approach  in  section  4.1.  The  final  process  to  "assign 
SNR  and  threshold"  to  the  peaks  uses  the  algorithms  outlined  in  section  3.4. 

The  algorithm  based  on  Figure  13  was  quantitatively  examined  by  observing  the 
number  of  peaks  over  15  dB  foimd  in  a  data  set.  This  level  is  well  above  the  noise 
and  we  would  expect  any  peaks  found  to  be  caused  by  target  signals.  The  number  of 
peaks  was  counted  for  the  pre-whitening  method  and  for  the  standard  post¬ 
whitening  method.  The  difference  between  the  two  (post-whitening  method 
subtract  pre-whitening  method)  is  plotted  in  Figure  14  below.  Positive  numbers 
indicate  that  the  pre-whitening  me&od  foimd  less  peaks  than  the  post-whitening 
method.  The  data  used  for  this  analysis  was  from  a  single  radar  region  over  Darwin 
airport.  The  region  contained  several  targets  as  well  as  the  Darwin  beacon. 
Propagation  was  present  in  two  ionospheric  layers  so  there  were  three  propagation 
modes  for  each  target.  This  produced  a  large  number  of  signals  above  15  dB. 


Figure  14  Improvement  by  peak  detection  on  pre-whitening  data 

As  can  be  seen  from  Figure  14,  the  pre-whitening  method  reduces  the  number  of 
multiple  peaks  on  average.  One  unexpected  result  was  that  the  pre-whitening 
method  gave  more  peaks  in  about  10%  of  cases.  These  were  not  investigated 
because  of  the  low  number  but  the  results  show  that  the  whitener  occasionally 
smoothes  the  data,  thereby  suppressing  peaks. 
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5.  Peak  Data  Structures 

The  final  issue  addressed  is  the  modification  of  the  peak  detection  scheme  to  use 
data  structures  that  incorporate  greater  flexibility.  The  current  peak  detection 
scheme  uses  fixed  array  dimensions  to  store  all  of  ^e  details  about  each  peak.  This 
structure  is  simple  and  suitable  for  an  array  processor  that  does  not  have  dynamic 
data  storage.  The  array  is  made  reasonably  large  to  cope  with  areas  of  high  peak 
density  and  has  space  allocated  for  N  peaks  in  each  range-azimuth  cell.  Overall,  the 
number  of  peaks  found  in  a  dwell  of  data  is  commonly  around  one  himdred.  This  is 
about  one  eighth  of  the  space  allocated  in  the  array. 

The  Stage  B  tracker  assumes  an  a  priori  distribution  for  false  alarms  with  a  preset 
detection  threshold  for  peaks.  The  Stage  B  peak  detection  scheme  described  in 
Section  2  has  an  unknown  higher  detection  threshold  when  the  number  of  peaks  per 
range-azimuth  cell  exceeds  the  upper  limit.  If  the  upper  limit  for  the  niunber  of 
peaks  stored  is  N  and  the  number  of  peaks  foimd  is  M,  the  above  situation  arises 
when  M>N .  In  this  case  we  have  (M-N)  peaks  above  the  threshold  that  are 
discarded.  This  is  equivalent  to  increasing  the  threshold  by  an  unknown  amoimt  for 
the  range-azimuth  cell  that  has  the  M  peaks. 

This  effect  can  be  beneficial  when  a  large  number  of  peaks  are  formed  due  to  meteor 
activity  or  spread  clutter  as  the  number  of  false  peaks  in  the  area  is  limited.  This 
prevents  these  areas  from  capturing  an  excessive  number  of  peaks  at  the  expense  of 
areas  with  less  clutter.  However,  it  does  invalidate  the  tracker's  a  priori  tlueshold 
assumption. 

Since  array  processors  are  no  longer  used  for  IRSU,  the  fixed  data  structures  are 
replaced  by  a  scheme  that  has  a  fixed  upper  size  for  the  number  of  peaks  in  each 
dwell.  This  allows  a  reduction  in  the  amount  of  data  transmitted  between  machines. 

The  fixed  upper  limit  is  achieved  by  adapting  the  detection  threshold  to  keep  a  near 
constant  number  of  peaks  in  each  dwell  with  a  variable  number  of  peaks  per  range- 
azimuth  cell.  This  leads  to  a  known  threshold  that  is  viniformly  applied  across  five 
whitening  ARD  data.  This  system  also  allows  the  implementation  of  other 
thresholding  techniques  such  as  a  higher  threshold  in  cluttered  areas.  What  is 
important  is  that  the  tracker  now  has  accurate  knowledge  of  the  threshold  for 
probability  calculations  imder  all  conditions. 

To  remove  the  problems  described  above  it  is  necessary  to  allow  the  number  of 
peaks  stored  in  a  particular  range  azimuth  cell  to  be  dynamically  altered.  Therefore, 
the  data  structure  used  for  peak  storage  needs  to  be  more  adaptive.  We  will  now 
consider  possible  new  data  structures. 

5.1  Linked  Lists 

We  could  optimise  the  space  used  in  peak  storage  by  using  a  linked  list  of  peaks  (see 
appendix  A).  In  this  case  we  would  be  allocating  only  enough  space  as  was  required 
by  the  peaks  that  were  foimd.  We  would  not  need  to  limit  the  number  of  peaks  to 
any  maximum  value  because  we  would  be  using  d5mamic  memory  allocation. 
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5.1 .1  Shortcomings  with  Linked  List  Approach 

The  linked  list  approach  is  not  appropriate  for  use  within  the  IRSU  signal  processing 
structure  for  several  reasons.  The  first  reason  is  an  incompatibility  with  Fortran. 
Fortran  is  the  language  used  for  the  peak  detector  code.  Fortran  does  not  have 
suitable  pointer  structures  and  d)mamic  memory  allocation.  This  is  not  a  major 
problem  since  we  can  write  the  list  functions  in  a  more  suitable  language  such  as  C 
and  call  the  functions  from  the  peak  detection  code,  however,  it  adds  to  the  overall 
completity  of  the  software. 

A  more  serious  problem  is  that  a  list  structiue  requires  a  set  of  pointers  into  the  list 
to  identify  which  peaks  belong  to  particular  range  azimuth  cells.  The  tracker 
requires  tirese  pointers  for  the  association  of  peaks  to  tracks.  This  requirement  is 
best  achieved  by  creating  an  array  of  lists  rather  than  a  single  list. 

The  third  problem,  however,  is  more  fundamental.  Linked  lists  are  pointer  based 
structures  and  as  such,  the  only  data  stored  in  the  array  would  be  pointers  to 
memory  locations  where  the  peaks  are  stored.  In  the  IRSU  signal  processing  suite, 
the  stage2  and  stage3  processes  run  on  different  machines.  This  means  that  if  we 
passed  the  linked  list  and  pointers  from  stage2  to  stageS,  the  receiving  task  would 
locate  this  data  in  different  memory  locations  in  the  stage3  processor.  Hence  the 
pointers  would  be  pointing  to  the  wrong  memory  locations  in  the  stage3  processor. 
To  avoid  such  problems,  we  would  need  to  "vmdo"  the  linked  list  by  packing  all  the 
peaks  together  and  then  sending  this  array  to  stage3.  In  stage3  we  would  need  to 
reconfigiu-e  the  linked  list  and  pointers.  This  is  all  a  bit  pointless  as  we  could  just 
use  the  packed  array  for  all  steps. 

These  difficulties  are  a  feature  of  the  pointer  nature  of  the  structure  and  so  all  other 
pointer  structures  (such  as  trees)  will  be  faced  with  ti\em. 

5.2  Packed  Array  Structure 

One  way  that  we  can  implement  a  list  structure  without  using  pointers  is  to  use  a 
fixed  size  array.  Each  new  data  point  is  entered  into  the  array  after  the  last.  We  do 
not  need  pointers  to  following  data  entries  because  we  can  find  the  next  entry  by 
incrementing  an  array  index.  This  structure  has  a  few  shortcomings  also. 

When  compared  with  a  linked  list,  deletion  of  data  from  the  array  becomes  a  very 
computationally  intensive  procedure.  If  data  is  deleted  from  the  middle  of  the  array, 
then  the  end  of  the  array  must  be  shuffled  along  to  fill  the  space.  This  is  required  to 
maintain  the  correct  ordering  of  the  list  elements.  In  a  linked  list  it  would  be  a 
simple  matter  of  reallocating  one  pointer.  Similarly,  adding  new  data  to  the  middle 
of  the  list  becomes  difficult  (the  same  shuffling  is  required). 

Because  we  are  not  using  the  array  as  a  true  data-base  we  can  ignore  the  problems 
above.  Data  will  never  need  to  be  deleted  from  the  array  (peaks  are  currently 
canceled  by  giving  them  cancel  values).  If  we  process  the  data  in  a  meaningful  order 
and  always  add  peaks  to  the  end  of  the  list,  the  peaks  will  be  stored  in  that  same 
order.  To  allow  for  a  variable  threshold,  a  temporary  working  array  is  used.  Once  a 
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smtable  threshold  has  been  decided,  the  peaks  above  this  threshold  are  moved  into 
the  output  peak  array. 

5.2.1  Indexing  the  array 

Like  the  linked  list  structure,  the  list  array  needs  some  additional  information  to 
allow  us  to  identify  peaks  in  some  given  area.  This  is  done  by  creating  an  array  of 
peak  indices.  In  the  stage3  processing,  the  tracker  requires  to  access  peaks  based  on 
the  range-azimuth  cell  in  which  they  are  detected.  To  enable  efficient  access  to  the 
peaks  the  peak  indices  array  stores  an  index  into  the  peak  list  array  where  the  first 
peak  for  a  particular  range-azimuth  cell  is  stored.  The  array  also  stores  the  number 
of  peaks  foimd  in  each  range-azimuth  cell.  If  no  peaks  are  foimd  in  a  particular  cell, 
both  of  these  fields  are  set  to  zero.  Combining  these  two  structures,  we  now  have 
the  capability  of  applying  a  variable  threshold  to  the  whole  dwell  whilst  not  limiting 
the  number  of  peaks  in  any  part  of  the  dwell. 

The  peak  list  array  and  the  peak  index  array  are  show  below  in  Figure  15. 
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PEAK  INDICES  ARRAY 
Figure  15  Peak  detection  data  structures 


Each  of  the  entries  in  the  peak  list  array  contains  information  about  a  peak.  This 
information  is  the  position  of  the  peak  in  azimuth,  range  and  Doppler  and  also 
additional  characteristics  such  as  peak  SNR. 
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6.  Implemented  Peak  Detection  Improvements 

The  changes  detailed  in  the  previous  sections  have  been  brought  together  to  produce 
the  enhanced  peak  detector.  The  enhanced  peak  detector  has  also  been  integrated 
with  the  adaptive  clutter  model  developed  by  Branko  Ristic  [3].  This  clutter  model 
develops  a  map  in  azimuth-range-doppler  space  that  divides  each  data  set  into 
clutter  zones.  The  probability  density  ftmctions  for  those  zones  are  recursively 
estimated  using  histograms  of  the  peaks  found  in  each  clutter  zone.  Part  of  this 
process  is  performed  during  peak  detection  and  part  of  it  is  performed  by  d\e 
tracker. 

When  a  new  dwell  of  data  arrives,  the  first  step  is  to  form  the  clutter  map.  This  is  a 
multidimensional  array  that  classifies  every  point  in  azimuth-range-doppler  space  as 
one  of  a  number  of  clutter  zones.  This  mapping  procedvue  uses  the  whitening  array 
as  a  starting  point  and  uses  various  image  processing  techniques  to  produce  the 
clutter  map  [3]. 

Once  the  clutter  map  is  formed,  peak  detection  is  carried  out  and  the  map  is  used  to 
associate  every  peak  foimd  with  one  of  the  clutter  zones.  Relevent  clutter  zone 
statistics  such  as  zone  ARD  volume  and  the  total  number  of  peaks  in  each  zone  are 
calculted.  The  map  statistics  (not  the  map  itself)  are  then  passed  on  with  the  peak 
data  and  the  tracker  performs  the  remaining  parts  of  the  clutter  modeling.  Namely 
the  recursive  formation  of  clutter  zone  probability  density  functions  based  on  the 
distribution  of  peaks  in  each  zone. 

The  peak  detector  also  uses  the  surface  clutter  edges  formed  by  the  data  whitener  as 
boundaries  for  detection.  In  the  whitening  process,  the  Ordered  Statistic  Whitener 
identifies  the  edges  of  the  surface  clutter  for  each  range  azimuth  cell  using  the  clutter 
suppression  algorithm.  This  can  give  a  different  maximum  and  minimtun  clutter 
edge  for  each  range  azimuth  cell.  These  clutter  edges  are  used  to  prevent  the  peak 
detector  from  looking  for  peaks  in  the  stuface  clutter.  Because  this  has  an  impact  on 
the  peak  distribution  in  doppler  space,  the  clutter  edges  calculated  are  stored  in  the 
peak  index  array. 

The  changes  described  in  the  previous  section  result  in  a  very  different  method  for 
peak  detection.  The  modified  algorithm  based  on  the  flow  diagram  in  Figure  13  is 
described  below  and  an  expanded  flow  diagram  is  shown  in  Figure  16. 

The  new  method  for  peak  detection  is  as  follows: 

1.  Form  the  clutter  map  using  the  array  of  whitening  factors  for  each  sample 
point. 

2.  Using  the  whitening  data  array,  locate  positions  where  the  whitening  data 
passes  above  a  minimum  threshold. 

3.  At  each  of  these  threshold  crossings  determine  if  that  cell  corresponds  to  a 
local  maximum  in  the  pre-whitening  data. 

4.  For  each  local  maximum  use  the  pre-whitening  data  to  interpolate  the 
peak  in  position  and  determine  the  SNR  using  the  whitening  data. 
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5.  Save  each  interpolated  maximum  into  a  temporary  peak  array  and  update 
a  ciunulative  distribution  of  the  SNRs  of  local  maxima  for  this  dwell. 

6.  Determine  the  threshold  required  to  obtain  a  desirable  number  of  peaks 
for  tracking  and  save  this  value  for  later  use. 

7.  For  each  maximum  that  crosses  this  threshold,  save  the  interpolated  peak 
parameters  and  the  associated  clutter  zone  in  the  peaks  array. 

8.  Update  the  array  of  indices  into  the  peaks  array  for  use  by  tracking 

As  shown  in  Figiure  16,  the  data  passed  to  the  tracking  system  are:  peak  index  array, 
peak  list  array,  SNR  threshold,  and  the  clutter  map  statistics. 

Pseudo  code  for  the  algorithm  is  set  out  below. 

For  each  beam 
For  each  range 

For  each  doppler 

If  whitening  data  cell  is  above  minimum  threshold 
If  pre-whitening  data  cell  is  a  local  maximum 

Interpolate  peak  using  pre-whitening  data 
Store  peak  data  in  potential  peaks 
update  distribution  table 

endfor  {doppler} 
endfor  {range} 
endfor  {beam} 

Using  distribution  table,  select  threshold  such  that  the  total  number  of  peaks  <  N 

For  each  peak  in  potential  peaks 
If  peak  SNR  is  above  threshold 

Store  peak  data  in  peaks  array 
Update  peakjndices 
endfor  {peak} 

The  algorithm  flow  diagram  is  illustrated  below: 
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7.  Test  Results 


We  will  now  undertake  a  comparison  of  the  performance  of  the  two  peak  detection 
algorithms.  Firstly,  let  us  consider  the  peak  location  accuracy. 

7.1  Peak  Location  Accuracy 

It  is  difficult  to  measure  the  accuracy  of  the  peak  interpolation  procedure  with 
experimental  data  because  we  do  not  know  where  the  target  peaks  truly  lie. 
However,  for  this  analysis,  we  have  used  the  artificially  injected  calibration  signal  in 
the  IRSU  system.  This  signal  is  injected  in  the  middle  of  Doppler  space  at  a  fixed 
position  in  range  and  azimuth  for  each  dwell.  It  is  actually  injected  in  between  the 
cell  sample  points.  We  have  therefore  taken  a  scatter  plot  of  the  peaks  generated 
near  the  injected  position  of  the  calibration  signal.  The  plot  was  constructed  using 
the  peaks  declared  over  the  space  of  188  different  dwells  of  clean  noise  data.  That  is, 
the  data  was  not  corrupted  by  meteors,  lightning,  spread  clutter  or  radio  frequency 
interference. 

The  scatter  plots  for  the  two  algorithms  are  shown  in  Figure  17.  They  show  all  of  the 
peaks  foimd  near  the  injected  calibration  signal  and  are  plotted  as  range  against 
azimuth.  The  Stage  B  peak  detector  is  shown  on  the  left  and  the  modified  peak 
detector  on  the  right. 


Figure  17  Peak  scatter  around  injected  calibration  signal 

There  is  an  obvious  improvement  in  the  consistency  of  the  interpolated  position. 
Also  the  peaks  from  noise  are  more  imiformly  distributed. 

As  a  measm-e  of  the  accuracy  of  the  positioning  process,  the  variances  of  the  above 
scatter  plots  were  determined.  This  was  done  by  collecting  the  peak  from  each  dwell 
with  the  highest  SNR  and  calculating  the  sample  variance  of  the  time  history  of  these 
peaks.  The  sample  variance  is  calculated  as  follows: 
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var 


-  1 

where  M  =  sample  mean. 

^  1=1 


Table  7  shows  the  sample  mean  and  variance  of  the  calibration  returns  calculated  as 
above.  The  measurements  are  in  vmits  of  cells  and  cells  squared  respectively. 


Stage  B  Algorithm 

Modified  Algorithm 

Range:  Mean 

2.52 

2.51 

Variance 

8.2  X  10-3 

9.6  X  10-5 

Azimuth:  Mean 

1.90 

1.82 

Variance 

4.2  X  10-3 

2.3  X  10-4 

Doppler:  Mean 

64.87 

64.95 

Variance 

2.1  X  10-2 

3.4  X  10-4 

Table  7  Calibration  signal  location 

The  position  variance  has  been  reduced  by  a  factor  of  85  in  range  and  by  a  factor  of 
18  in  azimuth.  This  corresponds  to  factors  of  more  than  9  and  4  respectively  in  units 
of  standard  deviation.  This  is  clearly  a  very  large  improvement  in  position 
consistency. 

7.2  Interpolation  Bias 

Next  we  consider  the  bias  of  each  of  the  algorithms.  For  this  purpose,  we  use  the 
same  data  set.  There  is  a  single  target,  but  it  does  not  cause  any  noticeable  impact  on 
the  results.  Again,  all  of  the  peaks  were  collated  for  the  188  dwells  and  then  they 
were  plotted  in  range-azimuth  space.  A  scatter  plot  of  the  Stage  B  peak  detector  is 
shown  in  Figure  18  (a)  and  the  modified  detector  in  Figure  18  G?). 

It  is  apparent  from  Figure  18  that  the  new  algorithm  shows  no  significant  biasing  in 
peak  interpolation.  In  contrast,  the  previous  method  shows  a  considerable  degree  of 
clustering  of  peaks  about  the  cell  points.  This  figure  demonstrates  the  visual  impact 
of  the  reduction  in  interpolation  bias  from  all  the  enhancements  that  were  discussed 
in  this  paper.  Most  of  the  bias  reduction  is  from  data  compression  and  amplitude 
interpolation.  Figure  18  gives  a  visual  interpretation  of  the  histogram  variance 
results  expressed  in  Table  6. 

7.3  Tracking 

The  ultimate  improvement  in  performance  should  be  observed  at  the  tracker  output. 
However,  the  new  peak  detector  with  frie  adaptive  threshold,  integrated  clutter  map 
and  modified  data  structures  cannot  be  integrated  with  the  Stage  B  Tracking  system. 
Effort  to  measure  the  enhanced  peak  detector's  performance  improvement  on 
tracking  is  not  considered  worthwhile.  Instead,  the  new  peak  detector  has  been 
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integrated  with  the  Multiple  Model  Tracker  [2].  A  test  of  the  fully  integrated  signal 
processing  and  tracking  algorithms  has  been  carried  out  but  will  not  be  discussed  in 
this  document. 
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Figure  18  Comparison  of  peak  bias:  (a)  Stage  B  system,  (b)  New  system 
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8.  Conclusion 


A  new  algorithm  for  detection  of  peaks  in  IRSU  data  has  been  developed.  The  new 
algorithm  has  been  shown  to  exhibit  no  significant  bias  in  peak  position  location  and 
the  peak  location  variance  has  been  dramatically  improved.  The  algorithm  has  been 
developed  with  new  data  storage  structures  to  enable  more  flexible  detection 
schemes  and  non-uniform  peak  density.  The  peak  detection  threshold  has  been 
made  adaptive  to  cope  with  changing  clutter  conditions. 

The  new  algorithm  incorporates  the  following  enhancements: 

•  Peak  interpolation  for  data  in  dB  in  range,  azimuth,  and  doppler  for  reduced 
peak  bias 

•  Increased  peak  detection  neighbourhood  for  reduced  peak  splitting 

•  Peak  location  and  interpolation  using  pre-whitening  data  for  improved 
position  accuracy  and  reduced  peak  splitting 

•  A  dynamic  peak  threshold  for  adaption  to  varying  clutter  conditions 

•  Flexible  data  structures  that  support  this  adaption 

The  algorithm  has  also  been  incorporated  with  an  adaptive  clutter  modeling 
algorithm  and  integrated  with  the  multiple  model  track  initiation  filter.  Together 
with  these  algorithms  it  provides  an  advanced  system  for  target  detection  and 
tracking. 
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Appendix  A 

Linked  Lists 

The  linked  list  is  a  simple  data  structure.  Let  us  consider  the  problem  at  a  somewhat 
abstract  level.  We  have  some  data  which  we  wish  to  store  in  an  ordered  manner. 
Let  us  then  define  the  data  type  node_type  as  follows: 

structure  { 

data_type  data; 

key_type  key; 

node_type*  next; 

node_type*  previous; 

}  node_type; 

where  node_type*  denotes  a  pointer  to  a  record  of  type  node_type  using  C-like 
s)mtax. 

This  record  now  contains  four  elements.  Firstly,  we  have  the  data  that  we  wish  to 
store.  Next  we  have  the  key  which  is  used  to  order  the  data.  Finally,  we  have  two 
pointers.  The  idea  is  that  each  list  entry  occupies  a  separate  space  in  memory.  We 
maintain,  at  all  times,  a  pointer  which  points  to  the  first  node  and  by  using  the 
pointers  stored  in  each  node,  we  can  move  forwards  or  backwards  through  the  list. 

This  concept  is  illustrated  in  Figure  19  below. 


Figure  19  Linked  List 

It  is  now  a  trivial  operation  to  insert  a  new  entry  into  the  list.  We  simply  find  the 
place  for  it  and  change  the  associated  pointers.  This  is  in  contrast  to  the  use  of  an 
ordered  array  structure,  for  example,  in  which  we  would  have  to  shuffle  along  the 
portion  of  the  array  after  the  insertion  point  to  make  room. 

The  same  is  true  also  for  list  entry  deletions. 
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